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Abstract 

Exact solutions of Einstein field equations invariant for a non-Abelian 
2-diniensional Lie algebra of Killing fields are described. Physical proper¬ 
ties of these gravitational fields are studied, their wave character is checked 
by making use of covariant criteria and the observable effects of such waves 
are outlined. The possibility of detection of these waves with modern de¬ 
tectors, spherical resonant antennas in particular, is sketched. 


Introduction 

Gravitational waves, that is a propagating warpage of space time generated from 
compact concentrations of energy, like neutron stars and black holes, have not 
yet been detected directly, although their indirect influence has been seen and 
measured with great accuracy. Presently there are, worldwide, many efforts to 
detect gravitational radiation, not only because a direct confirmation of their 
existence is interesting per se but also because new insights on the nature of 
gravity and of the Universe itself could be gained. For these reasons exact 
solutions of the Einstein field equations deserve special attention when they are 
of propagative nature. The need of taking into full account the nonlinearity 
of Einstein’s equations when studying the generation of gravitational waves 
from strong sources is generally recognized m Moreover, despite the great 
distance of the sources from Earth (where most of the experimental devices, 
laser interferometers and resonant antennas, are located) there are situations 
where the non linear effects cannot be neglected. This is the case when the 
source is a binary coalescence: indeed it has been shown in that a secondary 
wave, called the Christodoulou memory is generated via the non linearity of 
Einstein’s field equations. The memory seems to be too weak to be detected 

*The article is dedicated, in the occasion of his 60th birthday, to Raphael Sorkin who first 
introduced gravitational particles with spin different from 2. 
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from the present generation of interferometers ESI ( even if its frequency is in 
the optimal band for the LIGO/VIRGO interferometers) but of the same order 
as the linear effects related to the same source, thus stressing the relevance of 
the nonlinearity of the Einstein’s equations also (soon) from an experimental 
point of view. 

On the theoretical side, starting from the seventy’s new powerful mathemat¬ 
ical methods have been developed to deal with nonlinear evolution equations. 
For instance, a suitable generalization of the Inverse Scattering Transform al¬ 
lows to integrate [ 7 | Einstein field equations for a metric of the form 

g = f {z, t) (^dt'^ — dz^) + hii (z, t) dx^ + ft.22 {z, t) dy^ + 2hi2 {z, t) dxdy. 

Indeed, the corresponding vacuum Einstein field equations reduce essentially^ 
to 


(aH-iH,^) ,^+(aH-iH,J ,5 = 0, 

where H = \\hab\\ , ^ = {t + z) /v^, y = {t — z) /v^, a = ^|detH|. This 
is a non-linear differential equation whose form is typical of two-dimensional 
integrable systems and can be integrated [3 by using a suitable generalization 
of the Inverse Scattering Transform, yielding solitary waves solutions. 

A geometric inspection of the metric above shows that it is invariant under 
translations along the x, y-ax.es, i.e. it admits two Killing fields, dx and dy, clos¬ 
ing on an Abelian two-dimensional Lie algebra A^- Moreover, the distribution 
T), generated by dx and dy, is 2 -dimensional and the distribution orthogonal 
to T> is integrable and transversal to T>. 

Thus, it has been natural to consider m the general problem of character¬ 
izing all gravitational fields g admitting a Lie algebra Q of Killing fields such 
that: 

I the distribution T>, generated by vector fields of Q, is 2 -dimensional; 

II the distribution T>^, orthogonal to T> is integrable and transversal to T>. 

According to whether dim^ is 2 or 3 , two qualitatively different cases occur. 

A 2-dimensional Q, is either Abelian (A2) or non-Abelian {02)- A metric g 
satisfying I and II, with Q = A2 or ^2, will be called Q -integrable. 

The study of A2-integrable Einstein metrics goes back to Einstein and Rosen 
EOl, Rosen m, Kompaneyets [23, Geroch ini, Belinsky and Khalatnikov (Sj. 
Recent results can be found in cni 

The greater rigidity of C/2-iiitegrable metrics, for which some partial results 
can be found in Eauini, allows an exhaustive analysis. It will be shown that 
they are parameterized by solutions of a linear second order differential equation 
on the plane which, in its turn, depends linearly on the choice of a j-harmonic 
function (see later). Thus, this class of solutions has a bilinear structure and, 
as such, admits two superposition laws. 

^The function / can be obtained by quadratures in terms of the matrix H. 
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When dim ^ = 3, assumption II follows from I and the local structure of 
this class of Einstein metrics can be explicitly described. Some well known exact 
solutions innun], e..g. Schwarzschild, belong to this class. 

Besides the new local ^ 2 -integrable solutions, a procedure to construct new 
global solutions, suitable for all such ^-integrable metrics, will be also described. 

The paper is organized as follows. In section 1 gravitational fields invariant 
for a two dimensional Lie algebra are characterized. In section 2 the Einstein 
equations for such metrics are reduced, by using the symmetry, to the so called 
/r-deformed Laplace equation. Harmonic coordinates are also introduced. Sec¬ 
tion 3 is devoted to the analysis of the wave-like character of the solutions 
through the Zel’manov and the Pirani criterions. In section 4 the canonical and 
the Landau energy-momentum pseudo-tensors are introduced and a compari¬ 
son with the linearised theory is performed. In section 5 realistic sources for 
such gravitational fields are described. Eventually, section 6 is devoted to the 
analysis of the polarization of the waves. 

In the following, Kil (g) will denote the Lie algebra of all Killing fields of a 
metric g while Killing algebra will denote a sub-algebra of ICil (g). 

Moreover, an integral (two-dimensional) submanifold of V will be called a 
Killing leaf, and an integral (2-dimensional) submanifold of T)^ orthogonal leaf 


1 Geometrical aspects 


• Semiadapted coordinates. 

Let g be a metric on the space-time M (a connected smooth manifold) 
and Q 2 one of its Killing algebras whose generators X,Y satisfy [-Y, K] = 
sK, s = 0,1 

The Erobenius distribution T) generated by Q 2 is 2-dimensional and a chart 
{x^, x^, x^, x'^) exists such that 

Erom now on such a chart will be called semiadapted (to the Killing fields). 

• Invariant metrics 

It can be easily verified 0211311331 that in a semiadapted chart g has the 
form 


g = gijdx^dx^ +2 (li + smix"^) dx^dx^ — 2midx^dx'^ + 

{x'^Y ~ dx^dx^ + 

2 (^ — sAx'^) dx^dx"^ Xdx'^dx'^, z = I,2;j = 1,2 
with gij, mi, li, A, /i, v arbitrary functions of (x^,x^). 
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• Killing leaves. 

Condition II allows to construct semi-adapted charts, with new coordi¬ 
nates such that the fields ei = djdx, e^ = d/dy, belong to 

In such a chart, called from now on adapted, the components Z^’s and 
rrii’s vanish. 

We will call Killing leaf an integral (2-dimensional) submanifold of V and 
orthogonal leaf an integral (2-dimensional) submanifold of . Since 
is transversal to V, the restriction of g to any Killing leaf, S, is non¬ 
degenerate. Thus, (S', gig) is a homogeneous 2-dimensional Riemannian 
manifold. Then, the Gauss curvature K (S) of the Killing leaves is con¬ 
stant (depending on the leave). In the chart {p = a;^|g, q = a;^|g) one 
has 

g|g = — 2s]lq + dp^ + 2 (jl — s\q^ dpdq + Xdq^, 

where A, ]l,v, being the restrictions to S of A, g, u, are constants, and 
K (S) = As^ . 

1.1 Einstein metrics when g{Y,Y) ^ 0. 

In the considered class of metrics, vacuum Einstein equations, = 0, can be 
completely solved 03|. If the Killing field Y is not of light type, i.e. g{Y, Y) ^ 0, 
then in the adapted coordinates {x,y,p,q) the general solution is 

g = f{dx^ ± dg^) -I- /3^[(s^fc^(7^ — 2slq + m)dp^ + 2{l — skq)dpdq + kdq^] (1) 

where f = — A± 0^/2s‘^k, and f3{x,y) is a solution of the tortoise equation 

(3 + A\Ti\f3 — A\ = u{x,y), 

where is a constant and the function it is a solution either of Laplace or d’ 
Alembert equation, A±it = 0, A± = ± dyy, such that {dxU0 ± {dyU0 0 0. 

The constants k,l,m are constrained by km — P = ^l,k 0 0 for Lorentzian 
metrics or by km — = ±1, /c ^ 0 for Kleinian metrics (if / > 0). Ricci flat 

manifolds of Kleinian signature appear in the no boundary proposal of Hartle 
and Hawking in which the idea is suggested that the signature of the space- 
time metric may have changed in the early universe. Some other examples of 
Kleinian geometry in physics occur in the theory of heterotic N = 2 string (see 
m and 0) for which the target space is four dimensional 

1.1.1 Canonical form of metrics when g{Y,Y) ^ 0 

The gauge freedom of the above solution, allowed by the function u, can be 
locally eliminated by introducing the coordinates (it, v,p, q), the function v{x, y) 
being conjugate to u(x,y), i.e. A±v = 0 and Ux = Vy,Uy = T^x- In these 
coordinates the metric g takes the form (local ’’Birkhoff’s theorem”) 
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9 = 


2s'‘ki3 


{du^ ± dv'^) + 0^[{s^k'^q^ 


2slq + m)dp^ + 2(1 — skq)dpdq + kdq^] 


with (3 (u) a solution oi (3 + A\n\(3 — A\ = u. 


1.1.2 Normal form of metrics when g(Y,Y) ^ 0. 

In geographic coordinates (d, p) along Killing leaves one has 

g\s=(3^ [dd^+T(d) V], 


where 3F (d) is equal either to sinh^ or — cosh^ d, depending on the signature 
of the metric. Thus, in the normal coordinates, {r = 2s‘^k(3,T = v,d,p), the 
metric takes the form 


5 = El 




+ e2r^ [dd'^ + T (d) dp'^] 


( 2 ) 


where £i = ± 1 , £2 = if with a choice coherent with the required signature 2. 

The geometric reason for this form is that, when g{Y, T) 7 ^ 0, a third Killing 
field exists which together with X and Y constitute a basis of so(2,1). The larger 
symmetry implies that the geodesic equations describe a non-commutatively 
integrahle system |42|. and the corresponding geodesic flow projects on the 
geodesic flow of the metric restricted to the Killing leaves. The above local form 
does not allow, however, to treat properly the singularities appearing inevitably 
in global solutions. The metrics o, although they all are locally diffeomorphic to 
©, play a relevant role in the construction of new global solutions as described 
later. 


1.2 Einstein metrics when g{Y,Y) = 0. 

If the Killing field Y is of light type, then the general Lorentzian solution of 
vacuum Einstein equations, in the adapted coordinates (x,y,p,q), is given by 

g = 2f{dx^ + dy^) + p[(w (x, y) — 2sq)dp^ + 2dpdq], (3) 

where p = A^ + B with kl, S G R, $ is a non constant harmonic function of 
X and y, f = (V$)^ y/\p\/p, and w {x,y) is solution of the p-deformed Laplace 
equation: 


pAw + Vp ■ Xw = 0. 

Metrics © are Lorentzian if the orthogonal leaves are conformally Euclidean, 
i.e. the positive sign is chosen, and Kleinian if not. Only the Lorentzian case 
will be analyzed and these metrics will be called of ( 1 / 2 , 2 ) — isotropic type 
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In the particular case s = 1, / = 1/2 and /i = 1, the above metrics are 
locally diffeomorphic to a subclass of the vacuum Peres solutions EH. that for 
later purpose we rewrite in the form 

g = dx^ + dy'^ + 2dudv + 2{ip^^dx + (pydy)du, (4) 


where 


u = , V = qe ^ + Lp{x,y,u), 

with ip{x, y, u) a harmonic function of x and y arbitrarily dependent on u. 

In the case yt = const, the ^-deformed Laplace equation reduces to the 
Laplace equation; for ^ = 1, in the harmonic coordinates system (x,y,z,t): 

z = [( 2(7 — u> (a;, y)) + e^]/2 

t = [{2q-w{x,y))e~^-eP]/2, 

the Einstein metrics lEJ take ESI the particularly simple form , 

g = 2f{dx‘^ + dy'^) + dz"^ — dt^ + d [w) d (In \z — t|]) ■ (5) 

This shows that, when w is constant, the Einstein metrics given by Eq. 
are static and, under the further assumption $ = xy/2, they reduce to the 
Minkowski one. Moreover, when w is not constant, gravitational fields look 
like a disturbance propagating at light velocity along the z direction on the 
Killing leaves (integral two-dimensional submanifolds of T>). 

2 Physical properties 

The wave character of gravitational fields has been checked by using co¬ 
variant criteria. In the following we will shortly review the most important 
properties of these waves. 

2.1 Asymptotic flatness 

A first step toward a physical interpretation of metrics © is to characterize those 
which are spatially asymptotically flat. For the metrics ®, in the vacuum case, 
we will consider (spatially) asymptotic flat a metric approaching the Minkowski 
metric for -I- ^ oo . In terms of the functions /, y and w, this asymptotic 

flatness condition reads: 

x^-I-^ oo => f ^ const, const, w ^ cix + C 2 y + c^, 

where ci, C 2 and C 3 are arbitrary constants and the behavior of w can be easily 
recognized by looking at the Riemann tensor. In order for metrics m to be 
spatially asymptotically flat p must be constant m For p = 1, the equation 
for w reduces to a two-dimensional Laplace equation, then in the vacuum case 
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the only possible choice is w = const, because the Laplace equation does not 
have solutions tending to a constant value. 

Let’s consider the non-vacuum case. The simplest source for metrics © is 
dust with density p and velocity with an energy-momentum tensor = 
pu^u,, [ni. When is a light-like vector field, this tensor can describe the 
energy and momentum of null electromagnetic waves. This is not surprising; 
in fact Peres himself mi ESI indicated this as a possible source for his metrics 
(PP-waves) which, as we know, are diffeomorphic to a subclass of solutions 
found in ESI El El- With this stress-energy tensor we can depict realistic 
astrophysical sources as Gamma Ray Bursts or Cosmic Strings^; the symmetries 
of the vacuum solution is preserved. 

Being the time coordinate in the Killing leaves, the dust will be chosen to 
move parallel to the light-like Killing field Y, i.e., with velocity = 6^'^. The 
non vacuum Einstein equations with the energy-momentum tensor 

— P fip^vp 

are fulfilled, with / and p the same as in the vacuum case and w solution of 
the p-deformed Poisson equation 

pAw \7p ■ \7w = 2fp'^p, 

with c = 1, 8ttG = 1. From the asymptotic flatness condition with p = 1, the 
equation for w reduces to a two-dimensional Poisson equation 

Aw = p. 

It is well known that, if p goes to zero fast enough, it is possible to find non 
trivial everywhere regular solutions w tending to a constant value. The function 
/ satisfies the equation 

/A/-(V/)' = 0, 

and this implies that the function 

tp = In I/I 

is harmonic. Thus, in order to have everywhere regular spatially asymptotically 
flat solutions, / and p must be constant functions and the fluid density p must 
tend to zero fast enough. 

However, if we admit 5-like singularity in the {x,y) plane {i.e., string-like 
singularity, by taking into account the third spatial dimension), spatially asymp¬ 
totically flat vacuum solutions with / 7 ^ const and w const can exist. In this 
limiting case in which p (x, y) —> 6 (x, y), the energy-momentum tensor becomes 
the one usually employed to describe the gravitational effects of cosmic strings. 

■^Possible observations of Cosmic Strings have been reported recently |4U|. 
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As far as Gamma Ray Bursts sources are concerned, they are usually mod¬ 
elled with null-like hypersurfaces^. In several papers (i.e. iniEni) the study 
of exact solutions, describing impulsive gravitational waves, have been posed 
starting from a metric tensor of a PP-wave and extended to quadratic curva¬ 
ture gravity EI1IS2 and a rough estimation for the intensity of such waves is 
provided. Assuming typical parameters for Gamma Ray Bursts like an energy 
flux ^ 10“^ erg ■ cm~^s~^, a temporal extension St ~ 10s and an energy density 
Pq ~ being the distance from the source and E the total energy emit¬ 

ted, the maximal amplitude for a gravitational wave signal arriving on heart 
would be of the order 2'Kpf^G5t^ ~ 10“^®, far below the sensitivity of modern 
detectors. 


2.2 Zelmanov’s and Pirani’s criteria 


To check the propagative nature of gravitational helds described by metrics 
© several covariant criteria have been employed. In the general case when / 
and pL are not constant functions the Zel’manov criterion m is satisfied m 
Moreover, when / is a constant function, the only nonvanishing components of 
the Riemann tensor held reduce to 


R 


txzx — 


'^ixx 

2(z-t)2’ 


R 


txzy — 


ixy 

2{z-ty' 


D _ 

^ty.y - 


( 6 ) 


which, w{x, y) being a harmonic function, are all harmonic functions of x, y. As 
a consequence, the generalized Zel’manov criterion is still satished |1 hj . 

Besides the Zel’manov-Zakharov criterion, the Pirani algebraic criterion is 
also satished. In light-cone coordinates {u = [z — t)l^/2, v = [z + t)/y/2), 
where the metrics given by Eq. © read 

g = 2f{dx^ + dy"^) + 2dudv + dw din |u|, (7) 


the vector helds du and are both isotropic. Moreover, the only non vanishing 
components of the Riemann tensor are 

Ruiuj — —d^jW 

and this corresponds to a type-lSi Riemann tensor in the Petrov classihcation. 
From the natural interpretation of the Pirani criterion m it follows that the 
gravitational wave propagates along the null vector held du, that’s to say the 
gravitational wave © propagates along the z—axis with the light velocity c = I. 


2.3 The energy-momentum pseudo-tensors 

The exact gravitational wave 

g = dx^ + dy'^ + dz"^ — dt^ d (w) d (In \z — t|), (8) 

® Null-like frequently considered as models for astrophysical processes involving relativistic 
jets and sudden acceleration of huge quantity of mass 







given by Eq. for = 1, f = lj2 , has the physically interesting form of a 
perturbed Minkowski metric with h = dwd\n |z — t| . Moreover, besides being 
an exact solution of the Einstein equations, it is a solution of the linearized 
Einstein equations on a flat background too: 




To study its energy and polarization, the standard tools of the linearized theory, 
and in particular the canonical energy-momentum pseudo-tensor, can be used 


[13 Eg. 


With h = d (w) d {\n\z — t\) the Tq component of the canonical energy- 
momentum tensor vanishes since h has only one index in the plane transversal 
to the propagation direction because the components of the tensor h cannot be 
expressed in the transverse-traceless gauge. 

The non vanishing components of the 4—momentum density tensor = 
Tq , denoting the Landau-Lifshitz energy-momentum pseudo-tensor , are 
j>'^and p^. 


/ = 4 (t - z)-" [Cl)2 + C2iw,,y)^] + 4 (t - z)-^ C3V-[|V«;|' Vrc], 

pi = 0 , p^ = 0 , 

p3 = 4 (t - z)“^ [Cl (w 2,3,)^ -I- C2(w,a:y)^] + 4 (t - z)“^ C3V-[|Vw;|^ Vw]. 


where Ci are some positive numerical constants, V = [dx, dy) and the harmonic- 
ity condition for w has been used 113 ESj. 

The use of the Bel’s superenergy tensor |g 



where the symbol * denotes the volume dual, leads to the same result. Indeed, 
the only non vanishing independent components of the covariant Riemann tensor 


RaP'yS — 9apR 


Ri313 — —W,11 ; R 1323 — —^C,12 ; I?2323 — —W,22 ■ 


It follows that the density energy represented by the Bel’s scalar 




the C/“’s denoting the components of a time-like unit vector field, depends on 
the squares of w,ij. 

Thus, both the Landau-Lifshitz pseudo-tensor and the Bel superenergy ten¬ 
sor single out the same physical degrees of freedom. In particular, we can take 
the components htx and hty as fundamental degrees of freedom for the gravita¬ 
tional wave ©■ 

Since p° = P^, these waves move at light velocity, according with the result 
obtained by the Pirani criterion. 
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2.4 Polarization 


Even more controversial than for the energy and momentum, the definition of 
spin or polarization for a theory, such as general relativity, which is non-linear 
and possesses a much bigger invariance than just the Poincare one, deserves a 
careful analysis. 

It is well known that the concept of particle together with its degrees of 
freedom like the spin may be only introduced for linear theories (for example 
for the Yang-Mills theories, which are non linear, it is necessary to perform a 
perturbative expansion around the linearized theory). In these theories, when 
Poincare invariant, the particles are classified in terms of the eigenvalues of 
two Casimir operators of the Poincare group, and W'^ where are the 
translation generators and = \e)j,vpGP'^ is the Pauli-Ljubanski polariza¬ 
tion vector with Lorentz generators. Then, the total angular momentum 
J = L -\- S is defined in terms of the generators as J* = ■ The 

generators P^ and span the Poincare algebra, XSO{3, 1) 

( -^pcr] = ~ ’Hpcr^vp ~ Vijp^ptr T Vva^pp) 

\ [Mpiy, Pp] = ii'HupPp ~ VppPv) (9) 

I [Pu,PA = ^- 

Let us briefly recall a few details about the representation theory of this al¬ 
gebra. The Pauli-Ljubanski operator is a translational invariant Lorentz vector, 
that is [Pp, Wy] = 0, [M^^.Wp] = i{ri^pWp — rj^pW,,). In addition it satisfies 
the equation 

WpP^ = 0 . ( 10 ) 

The unitary (infinite-dimensional) representations of the Poincare group fall 
mainly into three different classes: 

• P^ = > 0 , W'^ = —m‘^s{s -\- 1), where s = 0 , 1 ,... denotes the spin. 

From Eq. cm we deduce that in the rest frame the zero component of 
the Pauli-Ljubanski vector vanishes and its space components are given by 
Wi = \eiojkP ^so that This representation is labelled 

by the mass m and the spin s. 

• P^ = 0, = 0. In this case W and P are linearly dependent 


Wp = XPp- 


the constant of proportionality is called helicity and it is equal to ±s . 
The time component of W is IT° = P • J, so that 


A = 


PJ 


which is the definition of helicity for massless particles like photons. 
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• p 2 _ ^r 2 _ where p is a continuous parameter. This type of rep¬ 

resentation, which describes particles with zero rest mass and an infinite 
number of polarization states labeled by p, does not seem to be realized 
in nature. 

Let us turn now to the gravitational fields represented by Eq. © . As it has 
been shown, they represent gravitational waves moving at the velocity of light, 
that is, in the would be quantised theory, particles with zero rest mass. Thus, if 
a classification in terms of Poincare group invariants could be performed, these 
waves would belong to the class of unitary (infinite-dimensional) representations 
of the Poincare group characterized by = 0, = 0. Recall that, in order 

for such a classification to be meaningful and have to be invariants of 
the theory. This is not the case for general relativity, unless we restrict to a 
subset of transformations selected for example by some physical criterion or by 
experimental constraints. For the solutions of the linearized vacuum Einstein 
equations the choice of the harmonic gauge does the job @21 ■ There, the residual 
gauge freedom corresponds to the sole Lorentz transformations. 

The polarization of gravitational fields represented by Eq. 0 can be es¬ 
timated by looking at the transformation properties of the two independent 
physical components of this metrics under a rotation in the plane (a:, y) or¬ 
thogonal to the propagation direction. This physical components, htx and hty, 
have only one index in the {x, y) plane orthogonal to the propagation direction 
du- Under the infinitesimal rotation TZ in the (x, y) —plane they transform as a 
vector. 

Applied to any vector (vi,V 2 ) the infinitesimal rotation generator TZ , has 
the effect 

TZvi = V2 , TZV2 = —fl, 

from which 

TZ^Vi = -Vi 1 = 1 , 2 , 

so that iTZ has the eigenvalues ±1. Thus, the components of that contribute 
to the energy correspond to spin-1 fields. The reason why it’s commonly believed 
that spin-1 do not exist is that, in treating with the linearized theory, solutions 
are implicitly assumed to be square integrable: this is not the case for solution 
like metrics 0 - 

These solutions are interesting for (at least) two reasons. First, they are 
asympotically flat (at least with a ^-like singularity) in the plane transversal 
to the propagation direction. Second, they are both solutions of the linearized 
theory and of the exact theory, so that the spin-1 result cannot be attributed 
to the first order approximation. 

2.5 Raychaudhuri equation 

A possible approach to the analysis of physical effects of gravitational waves is 
based on the Raychaudhuri equation. This approach has the important advan¬ 
tage to be covariant so that it is not needed anymore to care about the choice 
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of a coordinate system. The evolution of a beam (congruence) of non light¬ 
like curves (i.e., trajectories of test masses which one would like to observe) on 
curved space-time is constrained by the Raychaudhuri equation: 

n2 

VaO = -RcbA^A^ + - 2(7^ - — -f Vb (V^^)^, (11) 

‘^ab = V[a^b], aah = ^ {aAb) — 6 {gah^ AaAb)/“i (y aA)^A}, (12) 

Va = e = VbA\ A.A'^ = Tl 

where Rcb are the components of the Ricci tensor, denotes the vector field 
tangent to the curve, the signature of g is (-|-,—,—,—), the upper and the 
lower signs in Eq. m correspond respectively to the space-like and time-like 
case. The functions 0, aab, ^ab are called the expansion, shear and twist of the 
congruence. 

In our case it is easier to perform computations in the ’’Peres system of 
coordinates” where the metrics read: 

ds^ = dx'^ + dy^ + (p{x,y,u)du'^ + 2dudv, = 0 (13) 

and, as a timelike curve describing the motion, it is convenient to choose 


so that 


du'P 


3 > 


2^ - 2a^ = 


2{yy^) 


y-2)^ 


Zy-2) 


3 ’ 


0 = - 


5 (pu'pf , {dx^f + {dyvY 


\2{ip-2Y 


y-2r 


(14) 


in which, assuming that — 2 > 0, it is possible to isolate the physical effects of 
the spin-1 gravitational waves described by Eq. It is clear that, due to the 
non trivial dependence of 6 on the transverse coordinates, spin-1 gravitational 
waves have distinguishing features with respect to spin-2 gravitational waves. 
These effects are manifest, for example, in a distribution of test particles that, 
in the case of spin-1 waves will experience a permanent displacement that is, to 
say, a memory effects m 


3 Detection of gravitational waves 

The observable effect of a gravitational wave acting on two nearby test masses in 
free fall is mathematically described by the Jacobi geodesics deviation equation. 
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If the separation four-vector between two test masses is the equation can be 
written in the following form: 

= E}t,,ZPv;^v- (15) 

where D denotes covariant derivation, s is the proper time along the reference 
geodesic, u'' = dx'^/ds is the velocity four-vector, is the Riemann tensor. 

It’s commonly believed that reasonable sources of gravitational waves are so 
far from earth that it’s always possible to consider the weak field approximation 
of Einstein’s field equations. In this approximation a gravitational wave is de¬ 
scribed by the perturbation h^i, to the Minkowsky metric. Moreover, according 
with the standard textbook analysis m, it is possible to choose a TT-gauge 
and express the geodesics deviation equation (CSJ in the following form: 


dt^ 2 ^ dt^ 


(16) 


where the geodesic deviation four-vector Z^^ has been chosen to be 



where X denotes the position vector of one of two particles in the comoving 
reference frame of the other. The r.h.s. of equation itTCI) appears as an effective 
Newtonian force acting on test masses. An essential feature of these equation is 
the possibility of factorizing the time dependence (see cni and (EIHl later). 

Now let’s consider Eq. m for metrics 0 in weak field and small velocities 
{ds ^ dt) limits. We get: 


vi y-.ij vk 

ISV = > = 1.2,3, (IT) 

Equations oa and GO) show that both spin-1 and spin-2 waves are transversal 
to the propagation direction. Equations 03 are less trivial than Eq.s (Cni). They 
cannot be integrated to give a general solution. Moreover, the dependence on the 
time variable cannot be factorized. A reasonable expectation is a permanent 
deformation on the initial distribution of test masses, as in the case of the 
Christodoulu memory cn, a well known effect due to the non-linearity of the 
gravitational field. Similar effects involved in physical process generated by 
burst sources are called burst with memory (BWM) [^. 

More in general, Eq. (11511 . for metric (@1) and Z^ = ( 0, A j , reads: 


= -g^>^X^djdhdu^. (18) 

The above equation shows that one can have either attraction or repulsion 
according to the choice of the function ip (x, y, u) which is constrained, outside 
the matter source, only by the condition to be a harmonic function of x and y. 
For example, the choice tp = p{x,y)a{u), with p a harmonic function of x and 
y and a a decreasing function of u, will give repulsion. This is not surprising 
because it is known from QFT that spin-odd bosons generate repulsion between 
particles of the same charge, the charge, in this case, corresponding to the mass. 
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3.1 Experimental devices 

In this section the possibility of detection of spin-1 gravitational waves from the 
experimental point of view is considered. Nowadays there are two kind of ground 
based instruments able to investigate gravitational waves in the high frequency 
band laser interferometers (IFOs) and resonant mass detectors 

m In spite of the high complexity of any experimental device (whether IFOs or 
resonant antennas), the detection principle is very simple. It essentially consists 
in a measurement of the displacement of test masses described by 03) and lig. 
In the following we will review some basic features of the interaction mechanism 
between such instruments and gravitational waves. 

In the case of IFOs the test masses are suspended mirrors and the displace¬ 
ment induced by gravitational waves is measured by laser interferometry. At 
present time the most sensible IFOs are the two LIGO (with coherent antenna 
patterns) and VIRGO (with antenna pattern coherent with the one of GEO). 
Despite of their higher sensitivity, IFOs have non isotropic antenna patterns, 
i.e. their sensitivity strongly depends on the relative orientation of the incoming 
wave and the plane of the IFO’s arms. Moreover, a single IFOs cannot perform 
spin measurements. We will not consider them further. 

In the case of resonant antennas the detector is considered as an elastic body 
bathed by gravitational waves. The response of the detector can be studied by 
making use of the classical theory of elasticity [2B1; the generic infinitesimal 
mass element constituting the detector can be considered as a test mass and 
the relative displacement produced by gravitational waves is measured via the 
normal modes of oscillation. At present, cylindrical detectors (or Weber bars) 
are worldwide spread. They are generally three meters long and two tons heavy 
cylindrical objects made of aluminium. The most important resonant bars are 
those belonging to the IGEC network: ALLEGRO, AURIGA, EXPLORER, 
NAUTILUS, NIOBE Like an IFOs a single Weber bar cannot perform 
spin measurements and have non-isotropic sensitivity. One could imagine a 
combined use of several suitably oriented antennas to obtain some information 
about spin. At present time this is not feasible because they are oriented with 
coherent antenna patterns to reduce the false alarm probability. Therefore, from 
the point of view of spin measurement, the whole array is equivalent to a single 
bar. In this contest we will focus our attention on spherical detectors because, 
unlike other detectors, in principle they are able to determinate the polarization 
of any incoming gravitational wave (and in a wider sense, to distinguish between 
different metric theories of gravitation m) so they seem to be the most natural 
instruments to investigate the spin-1 gravitational waves. 

In order to describe the effect of spin-2 gravitational waves on an infinitesimal 
mass element of the detector located at position r = in a reference 

frame whose origin is at the center of the detector, it will be useful (also for 
later purpose) to denote by ^he force in Eq. llTCll : 
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Factorizing the radial and the angular dependence in /s=2(^)^) obtains a 
result depending on spherical harmonics with I = 2 only : 

fs=2{'^,t) = ( 20 ) 

' m 

Let be the displacement vector of an infinitesimal mass element, 

located at position r with respect to the center of mass of the initially unper¬ 
turbed solid with constant density p and Lame coefficients A and p. When a 
force / is acting on the body, the induced displacement vector ,t) 

field is solution of the following system of partial differential equations^: 

= + + J (21) 

In the following we will describe by external force / (pr ,t) the effects com¬ 
ing from the geodesic deviations®, whose explicit espression is given by ra or 
equivalently CT . A generic solution of Eq. m is expressed as linear combina¬ 
tion of eigenfunctions, with time dependence appearing in the coefficients am(t) 
only: 


m 

where is eigenfunction of the equation 


= (A-H/r)V(V • Vm)(22) 

and describes a free oscillation with frequency cOm- This model applies to any 
resonant mass detector once we impose the boundary conditions determined by 
the detector’s shape, i.e. to cylindrical antennas in three dimensions m or in 
one dimensional approximation EH, to spherical detectors [numEzi. 

The boundary conditions for sphere’s surface free from stress and strain are 
expressed by the following relation: 

cTijUj = 0 at r := | V| = R, (23) 


R being the radius of the sphere, rij the components of the outgoing surface 
normal vector, aij = XukkSij + 2puij the stress tensor and Uij = ^ {ui,j +Uj,i ) 
the strain tensor. Then, the time-dependent normal mode amplitude am(t) is 
solution of a driven harmonic oscillator equation: 


= - / u 

pNm J 

^Dissipation terms can be easily included in realistic cases 

®In the context of resonant antennas / (Y, i) generally contains two contributions: one 
from geodesic deviations and one describing the forces between the surface of the elastic body 
and other objects eventually matched with it (i.e. resonant transducers suitably tuned to the 
natural frequencies of oscillation). 


m(V) • / (V ,f)d®r, (24) 


^m(t} d.iy^Q7n(^) 
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with Nm normalization constants. All the interaction with external world (only 
gravitation in this case) can be included in the r.h.s. as additional terms of 
effective force acting on every normal mode. With some straightforward cal¬ 
culations Eq. 03 with boundary conditions 12211 can be solved giving the 
following solutions : 

+ T ~^— Ltfj{~r) + —V X Li/jiT"), (25) 

^irr l^div-free ^div-free 

with wave numbers ki^r = + 2/i, kdiv-free = where ‘pi/r ) and 

are scalar functions solutions of Helmotz’s equation, L is the angular mo¬ 
mentum operator L = x V, Cq, Ci and C 2 are constants, their numerical 
value depending on the boundary conditions in the specific case. To get 
regular solutions in r = 0 the scalar functions and must take the 

form: 


= ji{qr)Yim{'n), ipi'r) = ji{kr)Yim{'n), 

where ji is a spherical Bessel function and Yi^ifn) a spherical harmonic. 

Normal modes of oscillation can be divided in two family: Thoroidal modes 
(Co = C2 = 0) and Spheroidal modes (Ci = 0). The latter can be expressed in 
the form: 

= Ani{r)Yim{'n)Tt - iBni{r)lt X LYimifn), (26) 

with Ani{r) and Bni{r) combinations of Bessel functions. 

Let’s go back to the effective force (01) acting on every normal mode. If we 
chose for / (T^, t) the expression (1201, the only non-zero integrals on the r.h.s. 
of Eq. (1241) will be those containing the eigenfunctions ifnim with spherical 
harmonics Y 2 mi 0 , p): 

'uim = [ai(r)'n + l3i{r)RV]Yim{p ,p). 

Functions ai (r) and (3i (r) are combinations of spherical Bessel function of order 
2 and determinate the motion in radial and tangential direction respectively. 
According with this standard analysis, only few modes of the spheres, the five 
quadrupolar modes [1 = 2), will be coupled to gravity. Moreover, several studies 
|28| show that a finite number of resonators opportunely tuned to this mode’s 
frequency, will suffice to completely solve the problem of deconvolution of the 
signal revealed by spheres. 

Suppressing the index I we can write the r.h.s. of Eq. (1^ as: 

Fm[t) = f llm- f s=2d^r =]-hm{t)^ MR, (27) 

^sphere ^ 

the constant 7 depending on the physical properties of the elastic medium. 

In the case of spin-1 gravitational waves relation ( 1201 ) cannot be used. If we 
want to take into account the whole interaction of spin-1 waves with spherical 
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detectors, we will be compelled to solve the integral in the r.h.s. of Eq. iI23J) 
making use of geodesic deviation m and m whose r.h.s. will be denoted by 
Let us consider the most simple case choosing the I = 0 normal 
mode (1^ . For I = 0 thoroidal modes are absent while the non vanishing 
contribution to spheroidal modes comes from the term Ani{r)Yimi n) n . 
Even if it cannot be explicitly evaluated, there are no reasons for integral on 
the r.h.s. of (1^ to vanish 

F^it) = f lino ■ ^ 0 

Jsphere 

This computation, which heavily depends on the choice of the harmonic function 
w, does not lead to a general solution as in the standard spin-2 case, but stresses 
the interaction between an acoustic detector and a spin-1 gravitational wave. 
Moreover, since Fq is non vanishing, spherical modes with I = 0 are activated 
even in General Relativity too. Thus, the standard identification of normal 
mode’s index I with the spin component of the driving gravitational wave is not 
ensured. A further step in this direction would be to test the coupling between 
normal modes and spin components in more general cases. 


4 Conclusions 

It is still deep-seated the belief that spin-1 gravitational waves cannot be present, 
not only in General Relativity but in every metric theory of gravitation. It has 
been shown that this erroneous belief derives from the implicit assumption of 
considering only square integrable solutions of the linearized Einstein equations: 
indeed it is not true that it is always possible to reduce to TT-gauge and to 
remove all the non spin-2 components by a gauge transformation. 

Once we accept that gravitational waves may have spin-1 and may be emitted 
by reasonable sources, it becomes important to define the experimental condi¬ 
tions necessary to observe their spin. The first concrete possibility could be the 
use of spherical detectors.. Due to their resonant spectrum, described by Eq. 
(P|l . spherical detectors are the ideal instruments for studying the polarization. 
According to the standard theory, spherical devices are prepared to observe 
only spin-0 and spin-2 waves. In such an analysis the {x, ?/)—coordinates depen¬ 
dence is not taken into account. Clearly it will be difficult to detected spin-1 
gravitational waves with instruments having dimensions smaller than typical 
length scale of spatial variation of the waves. This does not implies that spin-1 
waves do not exist at all or that it is not possible to conceive new experimental 
apparatus capable to measure their spin. 
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